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ABSTRACT 

We combine published photometry for the nuclear star clusters and stellar spheroids 
of 51 low-mass, early-type galaxies in the Virgo cluster with empirical mass-to-light 
ratios, in order to complement previous studies that explore the dependence of nuclear 
star cluster (NSC) masses on various properties of their host galaxies. We confirm a 
roughly linear relationship between NSC mass and luminous host spheroid mass, albeit 
with considerable scatter (0.57 dex). In order to translate this to an A/nsc — c relation, 
we estimate velocity dispersions from the virial theorem, assuming all galaxies in our 
sample share a common dark matter fraction and are dynamically relaxed. We then 
find that Mnsc ~ cr 2 - 73 ± - 29 j with a slightly reduced scatter of 0.54 dex. 

This confirms recent results that the shape of the Mcmo — o relation is different 
for NSCs and super-massive black holes (SMBHs). We discuss this result in the context 
of the generalized idea of "central massive objects" (CMOs). 

In order to assess which physical parameters drive the observed nuclear cluster 
masses, we also carry out a joint multi-variate power-law fit to the data. In this, 
we allow the nuclear cluster mass to depend on spheroid mass and radius (and hence 
implicitly on velocity dispersion), as well as on the size of the globular cluster reservoir. 
When considered together, the dependences on Ms p h and i?s P h are roughly consistent 
with the virial theorem, and therefore MnSC c 2 - However, the only statistically 
significant correlation appears to be a simple linear scaling between NSC mass and 
luminous spheroid mass. 

We proceed to directly compare the derived NSC masses with predictions for two 
popular models for NSC formation, namely i) globular cluster infall due to dynamical 
friction, and ii) in-situ formation during the early phases of galaxy formation that 
is regulated via momentum feedback from stellar winds and/or supernovae. Neither 
model can directly predict the observations, and we discuss possible interpretations of 
our findings. 

Key words: galaxies: nuclei - galaxies: elliptical and lenticular, cD - galaxies: pho- 
tometry - methods: statistical 



1 INTRODUCTION 

The presence of a central massive object in the nuclei of 
galaxies appears to be a generic by-product of their for- 
mation and evolution. Over the last decade, it has become 
clear that in galaxies with massive stellar spheroids (i.e. 
massive ellipticals or early-type spirals with bulges more 
massive than > 10 10 Mq), the nucleus is nearly always oc- 
cupied by a super-massive black hole that outweighs by 
far the stellar mass within the central few parsecs. On 
the other hand, in most galaxies without massive stellar 

* E-mail: nleigh@rssd.esa.int 



spheroids, such as disk-dominated spir als and low-mass ellip - 
ticals, an SMBH is only rarely present l|Satvapal et al.[ [2009). 
Instead, the central mass concentration is dominated by a 
compact and massive nuclear star cluster. In between the 
two regimes, there appears to be a "transitional" spheroid 
mass range (8 < log(Ms p h[ M©]) < 10) over which the 
two CMO types can co-exist, and have comparable masses 
jGraham fc Spitler|[2009l : iNeumaver fc Walcherll2012r ). 

It is still unclear, however, whether these two in- 
carnations of CMOs form through different physics, or 
are merely different evolutionary stages of a common 
formation mechanism. This question has triggered a 
very active area of astrophysical research, with many 
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observation al studies of the statis ti cs a nd properties 



of N S Cs dCarollo. Stiavelli fc Mack] Il998l; iBoker et al 
20021. 12004; iBalcells et all 120031; iGraham fc Guzman 



20031 ; IWalcheret al l |2005|. I2OO6I; fcote et all |2006| ) and 



SMBHs (|Satvapal et al.l 120081 ; IGallo et alj l201Ch in low- 
mass spheroids, their co-existence ([Graham fc Driver] 



120071 ; ISeth et all 120081 ; iGonzalez Delgado et all |200S| ). 
as well as a numb er of proposed theories abo u t the 



growth of NSCs |Capuzzo-Dolcetta fc Miocch 



Agarwal & Milosavli evid 1201 ll ; lHartmann et al 



2008 



2011 



Antonini et alj I201S ) and the mutual feedback b e tween 



NSC and SMBH (iMcLaughlin. King fc Navakshinl 120061 : 



iNavakshin. Wilkinson fc Kindboolil ) 



I n 2006, three independent studies JWehner fc Harris! 
120061 ; iRossa et ai]|2006l ; IFerrarese et al.ll2006ah pointed out 
that NSC masses appear to correlate with the masses 
of their host spheroids, confirming an earlier result by 
IGraham fc Guz man ( 200jl) . Thes e scaling relations are simi- 
lar to those observed for SMBHs (IFerrarese fc Merritfeood ; 
iGebhardt et aklbOOrtlHaring fc Rixll2004l ), i.e. Mcmc cx a p . 
However, because of the observational difficulties to ob- 
tain reliable NSC masses, and the relatively small range of 
spheroid masses over which NSCs dominate over SMBHs, it 
is still debated whether both types of CMOs follow the same 
scaling relation (i.e. have the same value of /3), or whether 
they depend in different ways on the properties of their host 

spheroid. 

While IFerrarese et all (|2006al ) suggested that NSCs 
smoothly extend the SMBH relations to lo wer masse s , with 
/? ~ 4 in both cases, a recent study by IGraham] (|2012l ) 
concluded that the two types of CMOs depend differently 
on the total spheroid mass: NSCs follow a relation with 
(3 = 1.57 ± 0.24 while for SMBHs, 4 < /3 < 5. 

In order to add to the current discussion, we re- 
examine photometry of the stellar spheroids and NSCs 
in 51 early-type galaxies, obtained with the Hubble Space 
Telescope ACS Virgo Clu ster Survey , and published in 
IFerrarese et~"aH J2q06bl) a nd lCote et alj |2006l ). The analysis 
of lFerrarese et aT (2006a) compared the luminous mass of 29 
NSCs to the total mass of their host galaxies, the latter be- 
ing estimated from stellar velocity dispersion measurements. 
Here, we follow a somewhat different approach by comparing 
the luminous mass of both NSC and galaxy spheroids, de- 
rived via color-corrected mass-to-light ratios. This not only 
yields more reliable stellar masses, but also improves the 
statistics by increasing our sample size to 51 si nce we can 
use ne arly the entire sample of NSCs identified in lCote et alj 
(l2006h . 

This paper is structured as follows. In §[2l we summarize 
the data sources for our analysis and describe our method 
to convert the published photometry into estimates of the 
stellar mass of both NSCs and host spheroids. We present 
a quantitative analysis of the resulting Mnsc — Ms p h rela- 
tion in § [21 as well as of the inferred Mnsc — o relation. In 
addition, we demonstrate through a multi-variate analysis 
that the only statistically significant correlation with any 
observed parameter is a simple linear scaling of NSC mass 
with spheroid mass. In § [4] we briefly summarize the theo- 
retical predictions for two models of NSC formation, com- 
pare the present-day NSC masses to those predictions, and 
discuss the implications for the plausibility of either model. 
We conclude with a discussion and summary of our results 



and the implications for the origin of NSCs in § [5] and § [B] 
respectively. 



2 THE DATA 

For our analysis, we use data for 51 nucleated early-type 
galaxies observed during the Adv anced Camera for Surveys 
Virgo Cluster Survey (ACSVCS. ICote et al-lliool ). Taking 
advantage of the high spatial resolution of the Hubble Space 
Telescope, this data set provides detailed and accurate pho- 
tometry for both the galaxy bodies and their NSCs. We 
reject five galaxies from the original sample of ICote et al.l 
(2006) in which the apparent NSCs are significantly offset 
from the ga l axy's photocentre and therefore, as discussed by 
ICote et all |2006l ). may well be globular clusters that only 
appear to reside close to the nucleus due to a chance pro- 
jection. 

We note here that, given the morphological types of 
the sample galaxies (E, SO, dE, dSO, and dE,N), their stellar 
spheroids can be expected to be virialized, and to have little 
or no current star formation activity. This justifies use of a 
single (color-dependent) mass-to-light ratio for each galaxy 
spheroid in order to derive its stellar mass. 

We make use of the apparent z-band magni tudes, (g-z) 
color s, and half-light radii for b oth NSCs (fromlCote et al.l 
120061 ) and host spheroids (from IFerrarese et al. I l2006bl ). In 
order to convert to absolu te magnitudes and physical radii, 
we follow the approach of ICote et al.l |2006l ) and assume a 
common distance of 16.5 Mpc for all galaxies. For conve- 
nience, the observed galaxy and NSC magnitudes are listed 
again in Columns 2 and 3 of Table [1] We also list in Col- 
umn 4 the (completeness-corr ected) total n umbe r of globular 
clusters (GCs) as provided in lPeng et al] l|2008l . their Table 
2) which we will use in § 14.21 

In order to obtain estimates for the stellar masses of 
both NSCs and spheroids, we need to multiply their respec- 
tive z-band luminosities by an appropriate mass-to-light ra- 
tio. We use the empir i cally calibrated mass-to- light ratios 
provided bv lBell et alj |2003), accounting for the (g-z) color 
of NSC and spheroid, respectively, and list the resulting val- 
ues for Msph and Mnsc in Columns 5 and 8 of Table [T] 
We also provide estimates for the average GC mass ifiGC 
for our sample of galaxies. This is done by multiplying the 
z-band galaxy masses by the percentage of the total mass 
of the galaxy in GCs (obtained in the z-band) taken from 
IPeng et alj |2008l ). This gives an estimate for the total mass 
of the GC system, which we then divide by the total num- 
ber of GCs in order to arrive at an estimate for the average 
GC mass. Note t hat some of the percentages provided in 
IPeng et alj (|2008h are negative. In these few cases, we leave 
blank the corresponding entry for mac in Table [T] 



3 A/nsc SCALING RELATIONS 

In this section, we present a quantitative analysis of the 
observed Mnsc — Ms p h relation using estimates for lumi- 
nous masses derived from photometry. We also consider the 
inferred Mnsc — cr relation found by assuming virial equilib- 
rium, in addition to a multi-variate analysis that quantifies 
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Table 1. Properties of th e sample galaxies, iden tified by th eir VCS number (Column 1). Columns 2 and 3 list the absolute z-band 
magnitudes of the galaxy llFerrarese et alj|2006bl) and NSC l lCote et al.ll2006T ) (after accounting for a dis tance of 16 . 5 Mp c to get the 
absolute magnitude), respectively. Column 4 provides the completeness-corrected number of GCs (Nqq) jPeng et alj|200q ). Column 5 
provides the total stellar mass of the galax y Ms u h calculate d from its z-band magnitude; in units of 10 9 Mq. Column 6 gives the effective 
radius iisph of the galaxy spheroid {Ferrarese et al. 2006b). Column 7 lists the calculated velocity dispersion within t he effective radiu s 
of the galaxy in kms" 1 . Column 8 gives the inferred stellar mass of the NSC Af NSC , column 9 its half-light radius llCote et alj|2006l) . 
and column 10 its calculated velocity dispersion. Finally, Column 11 gives the average globular cluster mass rfiQQ in units of 10 s Mq. 
All listed errors denote the la uncertainties. 



VCS M z M z N GC M Sp h Rsph cr MnSC -%SC °"NC m GC 

(Galaxy) (NSC) (10 9 M ) (kpc) (kms" 1 ) (10 5 M ) (pc) (kms" 1 ) (10 5 M ) 



I 1 ) 


( 2 ) 


W 


W 


fx} 
(■>) 


/ p. \ 


(') 


(8) 


(9) 


(10) 




27 


-20.68 


-14.27 


62.0 ± 13.0 


23.59 


2.18 


241.19 


708.86 ± 53.55 


6.81 


236.61 ± 0.81 


4.57 


29 


-20.19 


-15.51 


84.0 ± 19.0 


16.08 


0.80 


328.87 


2188.37 ± 165.32 


25.87 


213.26 ± 1.03 


2.68 


30 


-20.73 


-13.49 


83.0 ± 25.0 


21.78 


2.00 


241.74 


251.25 ± 18.98 


1.92 


265.09 ± 28.65 


4.20 


31 


-20.29 


-12.95 


116.0 ± 24.0 


17.83 


1.34 


267.95 


231.43 ± 17.48 


2.80 


210.68 ± 13.11 


4.76 


32 


-20.26 


-14.41 


50.0 ± 14.0 


16.81 


1.09 


287.28 


862.06 ± 65.12 


12.89 


189.59 ± 3.04 


2.35 


36 


-20.0 


-14.13 


69.5 ± 9.8 


11.21 


0.84 


267.70 


623.10 ± 47.07 


16.66 


141.81 ± 2.97 


_ 


37 


-20.10 


-12.90 


20.1 ± 7.3 


9.88 


0.94 


237.25 


197.75 ± 14.94 


2.08 


225.96 ± 21.88 


2.95 


38 


-20.06 


-15.37 


47.0 ± 11.0 


13.93 


1.05 


267.09 


2118.26 ± 160.02 


40.13 


168.47 ± 1.03 


3.56 


39 


-19.64 


-17.14 


72.0 ± 12.0 


8.45 


1.29 


188.07 


8197.84 ± 673.67 


62.47 


277.04 ± 1.03 


4.58 


42 


-19.74 


-15.14 


71.0 ± 14.0 


8.74 


1.15 


201.91 


1615.18 ± 122.01 


47.81 


134.77 ± 4.30 


3.82 


44 


-19.44 


-11.58 


52.3 ± 8.5 


7.15 


0.89 


207.37 


38.14 ± 2.88 


4.81 


65.32 ± 1.34 


2.19 


47 


-19.70 


-12.02 


58.6 ± 9.3 


8.60 


1.56 


172.21 


89.90 ± 6.79 


4.24 


106.71 ± 3.02 


3.23 


48 


-19.21 


-12.81 


35.1 ± 7.6 


4.09 


1.61 


116.75 


142.51 ± 10.77 


2.88 


163.02 ± 9.69 


3.38 


49 


-19.71 


-12.48 


114.0 ± 12.0 


6.94 


1.22 


175.08 


113.26 ± 8.56 


3.28 


136.17 ± 6.48 


3.41 


50 


-18.73 


-12.09 


24.9 ± 6.0 


3.85 


1.76 


108.50 


73.43 ± 5.55 


2.80 


118.67 ± 7.38 


3.25 


51 


-18.57 


-10.79 


17.2 ± 5.4 


2.77 


1.01 


121.20 


21.05 ± 1.59 


2.48 


67.52 ± 5.07 


3.06 


52 


-18.32 


-13.61 


10.4 ± 5.0 


2.94 


0.55 


168.80 


304.45 ± 23.00 


8.65 


137.56 ± 0.74 


5.93 


55 


-18.33 


-12.46 


48.7 ± 8.4 


2.67 


0.96 


122.24 


100.97 ± 7.63 


3.04 


133.56 ± 7.26 


5.98 


57 


-18.25 


-13.24 


43.4 ± 7.9 


2.00 


1.34 


89.75 


196.63 ± 14.85 


12.25 


92.88 ± 1.38 


2.54 


58 


-18.27 


-9.90 


21.3 ± 6.1 


2.10 


0.74 


123.27 


7.76 ± 0.59 


1.76 


48.67 ± 5.90 


1.77 


59 


-18.16 


-11.02 


10.8 ± 5.6 


1.77 


2.44 


62.41 


23.81 ± 1.80 


3.04 


64.85 ± 3.52 


0.82 


60 


-18.64 


-12.20 


66.0 ± 9.5 


3.05 


1.63 


100.27 


88.16 ± 6.66 


2.16 


148.05 ± 13.60 


4.02 


62 


-18.34 


-12.03 


37.6 ± 7.4 


2.54 


1.46 


96.74 


60.80 ± 4.59 


9.53 


58.56 ± 0.49 


4.67 


63 


-18.19 


-9.89 


18.1 ± 5.5 


2.19 


1.46 


89.64 


10.74 ± 0.81 


15.70 


19.18 ± 0.38 


1.81 


64 


-18.25 


-12.55 


60.6 ± 9.3 


2.43 


0.79 


128.46 


104.93 ± 7.93 


18.66 


54.98 ± 1.25 


3.13 


65 


-18.04 


-9.78 


40.7 ± 7.6 


1.84 


0.79 


111.76 


7.21 ± 0.54 


1.44 


51.87 ± 8.13 


2.98 


66 


-18.32 


-9.86 


14.4 ± 5.7 


2.05 


2.17 


71.29 


10.60 ± 0.80 


1.92 


54.46 ± 5.88 


1.28 


68 


-18.69 


-12.08 


42.1 ± 7.9 


2.96 


2.08 


87.38 


61.80 ± 4.67 


6.65 


70.70 ± 0.31 


3.23 


69 


-18.18 


-11.89 


23.9 ± 6.1 


1.79 


1.53 


79.34 


56.71 ± 4.28 


2.32 


114.57 ± 9.50 


5.38 


70 


-17.01 


-9.82 


2.2 ± 4.2 


0.72 


0.73 


72.46 


7.21 ± 0.54 


2.56 


38.90 ± 2.79 


_ 


71 


-17.80 


-9.15 


15.5 ± 5.8 


1.10 


1.05 


75.13 


3.92 ± 0.30 


3.04 


26.31 ± 1.43 


0.36 


74 


-17.34 


-8.10 


7.3 ± 4.2 


0.70 


0.84 


66.84 


1.28 ± 0.10 


2.00 


18.57 ± 1.90 


0.29 


76 


-17.29 


-8.48 


6.3 ± 4.3 


0.87 


0.84 


74.74 


2.05 ± 0.16 


1.68 


25.62 ± 3.30 


0.96 


78 


-17.74 


-10.21 


54.2 ± 8.8 


1.33 


0.86 


91.24 


11.46 ± 0.87 


2.96 


45.59 ± 2.59 


2.51 


79 


-18.14 


-13.19 


-6.0 ± 13.0 


2.35 


0.57 


149.06 


197.79 ± 14.94 


9.69 


104.75 ± 0.93 




81 


-17.55 


-10.98 


16.5 ± 5.2 


1.11 


1.43 


64.56 


21.78 ± 1.65 


3.12 


61.23 ± 3.18 


2.21 


83 


-17.39 


-13.70 


3.6 ± 3.7 


1.19 


0.31 


143.46 


380.80 ± 28.77 


15.78 


113.91 ± 2.28 


39.88 


84 


-17.66 


-12.71 


26.7 ± 6.8 


1.31 


0.62 


107.05 


139.98 ± 10.57 


4.49 


129.54 ± 3.20 


1.43 


85 


-16.96 


-11.87 


28.7 ± 6.7 


0.52 


0.73 


61.78 


54.45 ± 4.11 


2.64 


105.24 ± 7.19 


9.82 


86 


-17.15 


-9.68 


9.2 ± 4.3 


0.78 


0.94 


66.97 


6.63 ± 0.50 


5.29 


25.96 ± 0.40 


1.79 


88 


-17.92 


-10.92 


9.2 ± 4.3 


1.79 


1.04 


96.26 


22.04 ± 1.66 


2.80 


65.01 ± 4.05 


2.53 


89 


-17.47 


-10.59 


20.4 ± 5.8 


0.97 


1.32 


62.89 


15.67 ± 1.18 


4.57 


42.96 ± 1.02 


4.52 


90 


-17.10 


-10.94 


0.6 ± 3.5 


0.36 


0.39 


70.05 


23.82 ± 1.80 


2.40 


73.00 ± 5.76 




91 


-17.43 


-11.69 


49.7 ± 8.6 


0.94 


0.87 


76.12 


42.52 ± 3.21 


10.17 


47.41 ± 0.48 


1.30 


92 


-16.93 


-10.05 


3.9 ± 2.6 


0.51 


1.07 


50.75 


9.60 ± 0.73 


3.28 


39.64 ± 1.89 


3.82 


93 


-16.94 


-12.70 


-9.0 ± 14.0 


0.80 


0.27 


127.22 


142.87 ± 10.79 


5.05 


123.38 ± 2.19 




95 


-17.12 


-11.28 


31.0 ± 7.0 


0.69 


1.79 


45.48 


32.10 ± 2.42 


21.22 


28.51 ± 0.70 


5.36 


96 


-17.37 


-11.18 


14.0 ± 5.7 


0.91 


1.69 


53.86 


26.00 ± 1.96 


4.00 


59.08 ± 1.90 


1.31 


97 


-16.71 


-12.18 


4.5 ± 3.9 


0.50 


0.55 


69.86 


78.02 ± 5.89 


2.00 


144.74 ± 14.80 


0.89 


99 


-16.41 


-9.58 


11.7 ± 4.8 


0.34 


0.89 


45.22 


5.61 ± 0.42 


4.65 


25.49 ± 0.58 


1.33 


100 


-17.40 


-11.84 


10.2 ± 4.7 


1.08 


3.15 


42.96 


51.42 ± 3.88 


6.57 


64.88 ± 0.32 


0.85 
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the dependence of Mnsc on the individual galaxy parame- 
ters. In all cases, the fits are normalized by dividing the pa- 
rameters being compared by their sample-averaged values. 
This significantly reduces the uncertainties in the offsets (i.e. 
the y- intercepts for our lines of best-fit). All uncertainties 
given in this paper correspond to the lcr-confidence level. 



3.1 The Mnsc — Ms p h correlation 

Figure [T] plots the stellar mass of the NSCs against that of 
their host spheroids. The solid line shows a weighted least- 
squares fit to the data found using a relation of the form 

Zo(?(Mnsc/10 7 ' 6 M ) = a + b ■ log(M Sph /W 9 - 7 M Q ). (1) 

The u ncertainties were calculated via a bootstrap method- 
ology (|Leigh. Sills fc Kniggei lioilh in which we generated 
1,000 fake data sets by randomly sampling (with replace- 
ment) data points from the observations. We obtained lines 
of best fit for each fake data set, fit a Gaussian to the subse- 
quent distribution and extracted its standard deviation. We 
find a = -0.53 ± 0.09 and b = 1.18 ± 0.16, with a scatter 
of 0.57 dex. A strong correlation between the two quanti- 
ties is evident, with a Spearman correlation coefficient of 
r s — 0.70, a significance level at which the null hypothe- 
sis of zero correlation is disproved of p s — 7.17 • 10 -9 , and a 
slope that is inconsistent with zero at the 7a confidence level. 
Table [T] does not list uncertainties for the spheroid masses 
provided in column 5, since we did not have uncertainties 
for the z-band magnitudes from which they were calculated. 
However, we performed a doubly-weighted least-squares fit 
using the spheroid m asses and uncertain ties provided in col- 
umn 8 of Table 1 in iPeng et al.l (|2008T ) (who derived their 
masses using different assumptions for the mass-to-light ra- 
tios than adopted in this paper). The slope and y-intercept 
for the resulting line of best-fit are the same as found using 
the spheroid masses calculated in this paper to within one 
standard deviation. 

This result is in broad agreement with various other 
studies. For example, using the HST imaging data of 
lLotz. Miller fc Ferguson! d2004T) for 45 nucleated d E,N g alax- 
ies in Leo, Virgo and Fornax. IWehner fc Harris] (2006) find 
values of a = -3.10 ± 0.81 and b = 1.18 ± 0.10. Us- 
ing HST near-infrared images of early-type spiral bulges, 
iBalcells. Graham fc Peletierl I (|2007f ) find a = -2.75 ± 0.15 
and b = 0.76 ± 0.13. Our analysis, which should be an im- 
provement over these earlier studies due to the larger sam- 
ple size and the use of color-dependent mass-to-light ratios, 
therefore re-affirms an approximately linear dependence of 
Mnsc on the luminous mass of the host spheroid. 

3.2 The Mnsc — cr correlation 

As mentioned earlier, we assume that the sample galaxies 
are in dynamical equilibrium, and that, therefore, the virial 
theorem can be used to predict the stellar velocity dispersion 
of the spheroid. W e use t he empirically calibrated formula 
of ICappellari et all (|2006h : 

^2 _ GMsph ^ 
5/ 9 J2s P h 

where A/s P h is the luminous mass of the galaxy spheroid, 
-Rsph is its half-mass radius, and f g = fi{,/f2 m is the bary- 



-0.5 

log (M. h /10« M 8 ) 



Figure 1. Correlation between the stellar mass of NSCs and their 
host spheroids. The solid line indicates the best-fitting linear rela- 
tionship of Equation[T] Err or bars were calcul ated using the 0.041 
mag uncertainty quoted bv lCote et all d2006T) ■ 



onic mass fraction. We assume that f g = 0.16 (|Spergel et al.l 
2003) for all sample galaxies, and that the measured half- 
light radii (Column 6 of Table [1} are a good approximation 
of the true half-mass radii, i.e. that luminous and dark mat- 
ter follow the same radial distribution. The resulting values 
for the predicted velocity dispersion are listed in Column 7 
of Table [T] Figure [2] then plots the relation between Mnsc 
and a. The best fit is again indicated by the solid line, it is 
described by 



log ( 



M- 



NSC 



10 7 - 6 M 



) = (-0.70±0.07)+(2.73±0.29)-/oc/(, 



128 km s 
(3) 

This correlation is slightly stronger (r s = 0.72 and p s — 
1.48 ■ 10~ 9 ) than that between Ms p h and Mnsc, and has a 
slightly reduced scatter (0.54 dex). 

Th e best-fit slope disagrees with that of iFerrarese et all 
(2006a) who conclude, based on measured velocity disper- 
sions of a smaller subset (29 objects) of the sample dis- 
cussed here, that log(M NSC ) = (6.91 ± 0.11) + (4. 27± 0.61) ■ 
Zog(o-/54kms -1 ). Taken at face value, this disagreement im- 
plies that our above assumptions are wrong, and that either 
the half-light radius is a poor indicator of the half-mass ra- 
dius, or that f g varies across the galaxy samp le. For example, 
in orde r for Equation[2] to yield the slope of IFerrarese et all 
(2006a), f g would have to scale with the square root of the 
luminous spheroid mass. 

On the other hand, it i s also possible t h at the relatively 
small sample used for the IFerrarese et alj l|2006al ) analysis 
causes systematic errors in the derived slope of the fit. This 
last poss i bility appears to be suppor ted by a recent s tudy o f 
iGraharr] (|2012l ) who re-analyzed the lFerrarese et all (|2006al ) 
data, together with complementary literature data, and con- 



!-)• 
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Figure 2. Correlation between the stellar mass of NSCs and the 
predicted velocity dispersion of their host spheroids. The solid line 
indicates the best-fitting linear relationship of Equation [3] Error 
bars are described as in Figure fT] 

eluded that Mnsc oc a 1 - 57±0 - 24 ^ m better agreement with 
Figured 

3.3 Multi-variate Analysis 

We now attempt to quantify whether or not the host galaxy 
properties are correlated in a way that is consistent with our 
use of the virial theorem as described in the previous section. 
More specifically, Equation[2] implies that the (total) mass 
of the spheroid and its half-mass radius are correlated. If the 
luminous matter indeed follows the same radial distribution 
as the dark matter, this should also be true for the luminous 
spheroid mass and its half-light radius i?sph- 

To assess whether or not this is indeed the case, we per- 
form a multi-variate analysis for our sample, using a gener- 
alized relation of the form 

M N sc ~ eM s a ph ^ ph iV GC , (4) 

where a, f3, 8, and e are all free parameters. 

Note that we have included the total number of GCs, 
Ngc, in this analysis, since the dependence of Mnsc on 
this parameter will be tested in §3] We perform a direct 
weighted least-squares fit to Equation [4] and calculate 2a 
and 3cr joint confidence intervals for the correlations be- 
tween possible pair of parameters a, /3, and 8. This is done 
by calculating the chi-squared value for every combination 
of each parameter pair, and identifying the minimum. We 
then find every combination that gives a chi-squared value 
that is within the specified confidence interval of the min- 
imum in order to obtain our joint confidence intervals. For 
a description of a similar application of this techn ique, we 
refer the reader to IVerbunt. Poolev fc Bassal \200^ ). 

Our goal is to determine which of the parameters in 



Figure 3. The large and small filled circles show the la and 
la, respectively, joint confidence intervals on the two parameters 
being considered for all three projections of the 4-dimensional pa- 
rameter space in Equation [4] The large solid circles and triangles 
show the exact power-law predictions for the GC infall and the 
in situ star formation models, respectively. 



Table 2. Best-fitting power-law indices for the parameters of 
Equation [4] 



Parameter Best-Fitting Value 

log e -i.o±i:g 

a 1 l +0 - 2 

" J - 1 -0.3 

-0.9t° i 



Equation [4] are correlated with Mnsc at a statistically sig- 
nificant level while naturally accounting for any correlations 
between them. In this way, we are able to quantify the de- 
gree to which Msph and ifeph are driving the correlation 
between Mnsc and a. 

The results are shown in Table [5] and illustrated in Fig- 
ure[3] As expected given the results of S I3.1I we find a roughly 
linear dependence of Mnsc on Ms p h- If the virial theorem 
is valid throughout our galaxy sample, Mnsc oc Ms p h im- 
plies that Mnsc oc R^K. The results for Rs p h appear to 
be consistent with this notion, but we caution that the un- 
certainty for this parameter is large. Therefore, our results 
suggest that Ms p h is the only galaxy parameter that corre- 
lates significantly with Mnsc- 



4 COMPARISON WITH PROPOSED 
FORMATION SCENARIOS 

In this section, we compare theoretical predictions for NSC 
masses from two popular formation models to the observa- 
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tional data. We begin by deriving theoretical relations for 
each formation scenario, which we use as a basis for com- 
parison with observational data. 



4.1 In-Situ Formation Regulated by Momentum 
Feedback 

A plausi ble mechanism for the birth of " s eed" N SCs, pro- 
posed by IMcLaughlin. King fc Navakshml 12006). is a n ex- 
tensio n of the arguments presented in lKind (|2003T ) and lKingl 
(2005) for the MsMBH-ff relation. It involves the formation 
of a massive star cluster from gas accumulating at the cen- 
ter during the early phases of the galaxy's evolution. Very 
soon after the onset of star formation, stellar winds from 
massive stars and supernovae drive a superwind from the 
nucleus. Initially, the wind is momentum-conserving with a 
momentum flux that is proportional to the Eddington lu- 
minosity. It sweeps up the surrounding medium into a thin 
supershell that gets pushed outward by the ram pressure 
it generates. Eventually, gravity balances the ram pressure 
and the shell stalls. More gas then falls into and feeds the 
nucleus. The process is repeated until the NSC reaches a 
critical mass, at which point the stall radius becomes very 
large (since the stall radius is proportional to the mass of 
the NSC). This gives the gas time to cool before it can reach 
the stall radius, since the cooling time begins to exceed 
the crossing time of the bubble. The gas becomes energy- 
conserving, and accelerates in its expansion to escape from 
the g alaxy. This fixes the mass of the NSC at the critical 
mass j|McLaughlin. King fc Navakshinll2006l ): 

Mksc=^, (5) 

where f g — 0.16 again is the baryonic mass fraction, k = 
0.398 cm 2 g _1 is the electron scattering opacity, and a is 
the velocity dispersion of the spheroid. If correct, this mech- 
anism would therefore n aturally explain the M nsc oc <t 4 
dep endence suggested by Ferrarese et al.l (|2006al ) . 

IMcLaughlin. King fc Navakshml (|200rj ) propose a value 
of A ~ 0.05 for the efficiency of massive-star feedback from 
a young nuclear cluster. This is based on simple calculations 
for the contributions from both supernovae and stellar winds 
to the total mo mentum flux, and agrees with more detailed 
calculations of iLeitherer. Robert fc Drissenl (|l992l ) for the 
total momentum deposition from solar- met allicity stars. 

We can re-write Equation [5] using the virial theorem. 
This assumes that the present-day mass distribution of the 
galaxy is representative of its distribution at the time of 
NSC formation. In particular, if the NSC formed before the 
galaxy mass was fully assembled, the velocity dispersion of 
the surrounding spheroid may well have been lower than it 
is today. 

Substituting the virial relation given by Equation[2]into 
Equation [5] this predicts for the NSC mass: 



Mk 



25/ 9 Att 

Re-writing this in a more convenient form 



M N sc = 1-66 x 10 b M Q ( 



(1\ / 2kpc \2/ M Sph \ 
VA/Viigph/ V1O 9 M / 



(6) 



(7) 



where we have assumed the baryonic mass fraction f g — 
0.16, and the electron scattering opacity k = 0.398 cm 2 g _1 . 



4.2 Infall and Merging of Globular Clusters 

The second NSC formation scenario occurs via the suc- 
cessive mergers of globular clusters (GCs) that spi- 
ral into the galactic centre due to dynamical friction 
(e.g. Tremaine. Ostriker fc Spi tzer 1975; Qu inlan fc Shapiro] 
1 199C ). Using N-body simulations of merging GCs in 
galaxy nuclei, several studies have successfully repro- 
duced some of the properties of NSCs, including their 
small effective radii and central v elocity dispersions (e.g. 
ICapuzzo-Dolcetta fc Miocchill2008l ). 

In general, the time required for a GC of mass wgc that 
starts on a circular orbit with radius r f rom the galactic cen- 
tre to spiral into the nucleus is given by ( Bi nnev fc Tremaind 
Il987l) : 



Tdyn [yr] = 



9 

To 



1.65 

G In A mGc 
3.73 x 10 11 
In A 



(8) 



V2kpc/ ^Okms-i/V m G c / 



v2kpc/ V250kms 



where a is the velocity dispersion of the galaxy. The 
Coulomb l ogarithm is denoted by In A, and is on the or- 
der of 10 (Bi nnev fc Tremaind 1 198*71 ). According to Equa- 
tion [8] dynamical friction operates most efficiently on the 
most massive GCs. 

If the NSC was indeed built from GCs that have sunk 
into the nucleus, its mass should follow a relation of the 
form: 



My 



/gcMgc ~ /GcJVGCfGC, 



(9) 



where Mgc is the total mass of the globular cluster system, 
Agc is the total number of GCs, mcc is the present-day 
average GC mass of the galaxy (Column 11 of Table[T]), and 
/gc is the mass fraction of the GC system that has merged 
onto the NSC. 

This mass fraction should be proportional to the ratio 
of the age of the galaxy, r ga i, to the time for an average GC 
to fall into the core due to dynamical friction, Tdyn- 



Mis 



/tot- 5 ^- Agc m G c. 

Tdyn 



(10) 



Here, the additional constant / to t accounts for the fact that 
only GCs that are located within a radius at which Td yn ^ 
Tgai will have had enough time to have fallen into the NSC 
via dynamical friction. 

Substituting Td yn = T ga i = 10 Gyr in Equation|8] shows 
that for most galaxies in our sample, only GCs well inside the 
half-light radius i?sph can be expected to have fallen into the 
nucleus. However, the present-day orbits of most GCs in our 
sample are compar able to or exceed the half-light radius in 
their host galaxies (|Peng et al.ll2008l ). In other words, most 
GCs will not have had sufficient time to have s piraled into 
the nu cleus, i.e. fee 1- The compilation of iPeng et al.l 
(2008) can be used to estimate the present-day value of /tot 
by calculating the fraction of GCs inside i?sph- This yields 
/tot ~ 0.2, in reasonable agreement with the notion that 
Agc has not changed significantly over the lifetime of the 
galaxy. 

We assume a universal initial GC MF that adheres to 
the same radial dependence in all galaxies, and use the dy- 
namical friction timescale calculated at the half-light radius 
as a proxy for the corresponding volume-integrated value. 
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Using Equation [8] we get from Equation 1 101 

M N SC = T S^ rfl GcRsph a ~ 1 / tot ' C 11 ) 

where we have used the half-light radius of the galaxy 7?s p h 
as a proxy for the average initial orbital radius from which 
clusters begin spiraling into the nucleus. 

Using the virial theorem to substitute for a (Equa- 
tion |5J, and /tot = 0.2 then yields 

Mnsc = 0.43(G/ 9 ) 1/2 hiAT sa im 2 GC M-^ 2 R-X 2 N GC . (12) 

Alternatively, we can re- write Equation [T^] in the more con- 
venient form: 

M NSC = 5.36-10 M ( — ) ( ){ W ^) 

(13) 

where we have assumed a galaxy age of 10 Gyrs, f g = 0.16, 
and In A = 10. 



4.3 Comparison with Observations 

We proceed to compare the predictions of the two formation 
scenarios discussed in the previous sections to the observed 
NSC masses, M£| s c . 

For the in-situ star formation model, the predicted NSC 
masses, M^q, as calculated from Equation [SJ| are plot- 
ted against the observations in Figure [4] A weighted least- 
squares fit yields the relation 

^(i<Sfc) = (0.68 ± 0.07) / O . g (Jg|-)- ( 0.24 ± 0.0 9 ), 

(14) 

The correlation is rather strong (r s = 0.75, p s = 2.34-10 10 ), 
and has a scatter of 0.54 dex. However, the offset be- 
tween the observed and predicted values suggests that the 
in-situ formation model over-predicts the observations by 
roughly two orders of magnitude across the galaxy sam- 
ple. Given that secular evolution process es such as addi- 
tiona l infa ll of cold molecular gas (e.g. ISchinnerer et all 
20061. 120071) and/or young star clusters may well have added 
to the observed present-day NSC masses, and that all mech- 
anisms for mass loss, such as the ejection of stars due 
to dynamical encounters, appear n egligible in NSCs (e.g. 
Spitzerl 1 19871 ; iHeggie fe Hutl 120031 ; IValtonen fe Karttunenl 
2006T) this appears to present a seve re challenge for the 
McLaughlin. King fe Navakshinl (|2006h model. 

Figure [5] makes the corresponding comparison for the 
GC infall model. Here, the predicted NSC masses were cal- 
culated using Equation 1121 The solid line in Figure [5] is the 
best fit linear relation between the two quantities, described 
by 

^iwfc) - (0.42 ± 0.26)/ O ,(^||-)-(0.74 ± 0.13). 

(15) 

The best-fit correlation has a slope that is marginally con- 
sistent with zero, has a large scatter of 0.80 dex, and is sta- 
tistically weak (r s = 0.49, p s = 5.24 • 10" 4 ). Moreover, the 
y-intercept suggests that the GC infall model under-predicts 
the observed NSC masses by about an order of magnitude 
across the galaxy sample. Adding to this problem is th e 
fact that, as pointed out bv lAgarwal fc Milosavlievid (|201l[ ). 
many GCs may be disrupted before reaching the nucleus, 
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Figure 4. Observed versus predicted NSC masses for the in situ 
star formation model, with the latter calculated using Equation[6] 
The solid line shows a weighted least-squares fit to the data ac- 
cording to Equation ll4l Error bars are described as in Figure [T] 



which would result in an even smaller NSC mass. We con- 
clude that the GC infall model provides a poor description 
of the observations. 

The multi-variate analysis presented in § 13.31 offers an 
additional test of the in situ star formation and GC infall 
models. In particular, from Equation [6] the in situ star for- 
mation model predicts a = 2, [3 = —2, and 8 = 0, whereas, 
from Equation 1121 the GC infall model predicts a = —1/2, 
P = —3/2, and 5 = 1. Therefore, if the various combinations 
of these values for a, /?, and 5 for one of our competing 
formation scenarios are found to lie within the joint con- 
fidence intervals in every 2-dimensional projection of this 
3-dimensional parameter space, this could be interpreted as 
evidence that the data are consistent with that particular 
formation scenario. 

The la and 2a joint confidence intervals for all three 
projections of our 3-dimensional parameter space are shown 
in Figure [3] The theoretically predicted relations are shown 
by solid triangles and circles for the in situ star formation 
and GC infall models, respectively, as described in the figure 
caption. The precise values for the power-law indices corre- 
sponding to both the in situ and GC infall models do not 
coincide with even the 2a confidence intervals in any of the 
three projections. 



5 DISCUSSION 

Our results suggest that Mnsc scales linearly with lumi- 
nous spheroid mass Ms p h- We also find a dependence of the 
form Mnsc oc a 2 - 73±0 - 29 when using velocity dispersions es- 
timated via t he virial theorem . This appears to support the 
conclusion of iGrahaml |2012i ) that the observed Mnsc — a 
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Figure 5. Observed versus predicted NSC masses for the GC 
infall model, with the latter calculated using Equation 1121 The 
solid line shows a weighted least-squares fit to the data according 
to Eouation ll5l Error bars are described as in Figure [T] 



relation is not an extension of the observed Msmbh — cr re- 
lation. Our results thus add further evidence to the notion 
that NSCs and SMBHs may not share a common origin, 
and that, instead, different formation mechanisms may be 
responsible for the two incarnations of a central massive ob- 
ject. We note that we find a slightly stronger correlation for 
the a — Mnsc relation than we do for the Ms p h — Mnsc 
relation. This raises the question whether perhaps a more 
fundamental relation is that between the velocity dispersion 
of the host spheroid and that of the NSC. 

Such a relation could be explained if the NSCs are in 
thermal equilibrium with their host spheroids, possibly be- 
cause galaxies tend to evolve towards (or be born in) such 
a state. To test this, we again used the virial theorem to 
calculate velocity dispersions also for all NSCs in the galaxy 
sampl e, using their effective radii measured by ICote et akl 
(2006), and listed in Column 9 of Tabled] Figure© shows 
that both quantities are indeed correlated, with a linear fit 
yielding 

log( T N ° 1 ) = (1.27±0.09)Zo 5 ( T ) -(0.01±0.02), 

y Vl02kms- 1 / v ; y Vl28kms- 1 / k ; ' 

(16) 

The correlation is rather strong (r s = 0.73), and has a 
smaller scatter (0.24 dex) around the best linear fit than 
either the Ms p h — Mnsc or the a — Mnsc relation. This 
supports the idea that the the dynamical evolution of NSCs 
is significantly affect ed by heat input from the surrounding 
galaxy (e.g. iMerritiJ 12009 ). If indeed the present-day NSCs 
have reached an energy balance with their host spheroids, 
it implies that the time-scale for energy exchange between 
them is short compared to the age of the galaxy. This sug- 
gests that the present-day properties of NSCs do not neces- 
sarily reflect those at the time of their formation, since there 



Figure 6. Comparison of the stellar velocity dispersion for NSC 
and host spheroid. The solid line shows a weighted least-squares 
fit to the data according to Equation ll6l Error bars were cal- 
culated us ing the 0.041 mag and 0.007 arcsecond uncertainties 
quoted by I Cote et al.l l|2006h for the NSC apparent magnitudes 
and radii, respectively. 

was ample time for them to be modified by the presence of 
the surrounding spheroid. 

The masses we have used for our analysis are derived 
from photometry, and therefore measure only the stellar 
mass of the galaxy. Measured velocity dispersions, on the 
other hand, provide a better proxy for the total mass of the 
spheroid, including the dark matter halo. It would therefore 
be most useful to perform a similar analysis using measured 
velocity dispersions once they become available. Neverthe- 
less, because of the arguments presented in Section \3. 31 our 
results have cast some doubt on the assumption that the 
baryonic mass fraction f g is the same for all early-type galax- 
ies and/or that the dark matter has the same radial profile 
as the luminous matter. 

We have also compared the observed present-day NSC 
masses to theoretical predictions from two proposed models 
for NSC formation, namely in-situ star formation regulated 
by momentum feedback and the successive infall of star clus- 
ters due to dynamical friction. As discussed in § 14.31 both 
models appear to have difficulties in predicting the observa- 
tions. 

The in-situ model predicts NSC masses that exceed 
the observed ones by about two orders of magnitude. Since 
the NSC may well have grown through secular processes 
throughout the life of the galaxy, this is probably only a 
lower limit for the true difference. The discrepancy may be 
somewhat reduced if one assumes that the host galaxy was 
significantly less massive and/or more extended at the time 
of NSC formation than it is today. In this case, Equation [6] 
would yield lower NSC masses. For example, if the NSC host 
was half as massive and twice as extended when the NSC 
formed, the resulting NSC mass should be ~ 16 times lower. 
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On the other hand, in this scenario it is difficult to explain 
why the Mnsc — cr re lation does not seem to have the pre- 
dicted slope of P = 4 l|Graharall20l3 y 

The GC infall model, on the other hand, systematically 
under-predicts the observed NSC masses by about an order 
of magnitude across the galaxy sample. This could, at least 
partly, be rectified by adopting a larger value for the average 
GC mass. Indeed, the most massive GCs are most responsive 
to dynamical friction, and therefore are most likely to have 
spiraled into the nucleus by today. Using the present-day 
average GC mass as an estimate for the average mass of 
the GCs that have merged onto the NSC in Equation [12] 
may therefore provide too low an estimate. Indeed, most 
NSC masses can be reproduced if 10 very massive GCs 
(> 10 6 M ) had fallen in via dynamical friction. With the 
data at hand, we cannot rule out this possibility. 

In summary, both models appear to have difficulties in 
explaining the observed NSC masses. Of course, this does 
not imply that neither in-situ star formation nor GC infall 
has occurred, but it suggests that either model by itself pro- 
vides an over-simplified description of NSC formation. 



6 SUMMARY 

In this paper, we have performed a statistical analysis to 
study the observed Mnsc — Ms p i> and inferred Mnsc — a 
relations using a large sample of 51 nucleated early- type 
galaxies taken from the ACSVCS. We used luminous galaxy 
masses, and calculated velocity dispersions found by assum- 
ing virial equilibrium applies for our sample of galaxies. We 
confirm a linear relationship between NSC mass and lumi- 
nous host spheroid mass reported by previous studies, but 
infer a s hallower slope o f the M nsc — a relation than re- 
ported in lFerrarese et all (|2006al ). 

Our results support a recent claim (Graham 2012) that 
NSCs do not obey the same M — a relation as do SMBHs, 
adhering instead to a shallower dependence. Specifically, we 
found Mnsc oc a 2 73 ± - 29 . This appears to cast doubt on the 
notion that NSCs and SMBHs share a common evolution- 
ary path, i.e. that they are merely different incarnations of 
"central massive object". We specu late that the dis c repanc y 
between our results and those of iFerrarese et all (|2006al L 
who found an Mnsc — o relation closer to A/nsc oc a 4 us- 
ing observed velocity dispersions, can perhaps be explained 
if the baryonic mass fraction is not the same for all early- 
type galaxies, or follows a different radial profile than the 
luminous matter. 

Finally, we have compared the observed present-day 
NSC masses to first-order predictions of two proposed 
NSC formation scenarios, namely GC infall and feedback- 
regulated in-situ formation. We find that neither scenario 
can provide a reasonable agreement with the observations. 
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